#ifndef __TENSOR_H
#define __TENSOR_H
#include "ELAPLA.h"
#define kroneker(i, j) ((i == j) ? 1 : 0)


double CalcuP(Stress sigma);

double CalcuQ(Stress sigma);

//Calculate Deviatoric Stress
//sigma:(in) Present stress
//Deviatoric: (out) Deviatoric stress
void CalcuDeviatoric(Stress sigma, Stress Deviatoric);


double J2(Stress sigma);


#endif